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Abstract. We introduce a new class of deterministic networks by as- 
sociating networks with Diophantine equations, thus relating network 
topology to algebraic properties. The network is formed by representing 
integers as vertices and by drawing cliques between M vertices every 
time that M distinct integers satisfy the equation. We analyse the net- 
work generated by the Pythagorean equation x 2 + y 2 = z 2 showing that 
its degree distribution is well approximated by a power law with expo- 
nential cut-off. We also show that the properties of this network differ 
considerably from the features of scale-free networks generated through 
preferential attachment. Remarkably we also recover a power law for the 
clustering coefficient. 

We then study the network associated with the equation x 2 + y 2 = z 
showing that the degree distribution is consistent with a power-law for 
several decades of values of k and that, after having reached a minimum, 
the distribution begins rising again. The power law exponent, in this 
case, is given by 7 ~ 4.5 We then analyse clustering and ageing and 
compare our results to the ones obtained in the Pythagorean case. 



1 Introduction 

The study of complex networks has recently attracted much attention in the 
physics community and currently represents an important area of multi-disciplinary 
research [T] . Network theory studies the properties of real world systems and phe- 
nomena displaying a web-like structure, usually represented mathematically as 
graphs. The individuals of the system under consideration are then symbolised 
as vertices and their relations or interactions as edges. For instance the World 
Wide Web is a network of web pages connected by hyperlinks; the cell can be 
described as a network of chemicals connected by chemical reactions; human 
language can be viewed as a web where words are linked if they appear adjacent 
or one word apart in sentences [112] . 

We recall that a network G is defined formally as a couple (V,E), where V 
is the set of vertices and E is a set of couples of vertices, called edges. If E is an 
ordered set, the network is said to be "directed". 

A network can be specified by its adjacency matrix A = {a^ }, where = 1 
if vertex i is connected to vertex j and Ojj = otherwise. The degree ki of 
vertex i is then defined as the number of nearest neighbours of vertex i, that 
is: ki = J2j a ij- The neighborhood Qi of a vertex i is the set of all vertices 
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connected to i. Notice that, in contrast with the terminology usually introduced 
in analysis and general topology, fli does not contain i as one of its elements. 
The degree function k% then "counts" the number of vertices belonging to each 
neighborhood and gives thus a first characterization of a network's topology. In 
many cases a more complete characterization of network topology is given by 
the weighted adjacency matrix W = {uiy}, where the numbers w,j represent 
the weight of the connection between vertex i and vertex j. This number can 
be defined in many ways. In our case we will define it as the total number of 
edges connecting vertex i and vertex j. We then define the strength Sj of vertex 
i as the total number of edges belonging to vertex i, that is Sj = J^ . Wij. When 
multiple edges are prevalent a weighted analysis is obviously worthwhile. More 
refined characterizations of network topology also take account of clustering p] 
and of the average distance between two vertices [1] . 

Complex networks were originally studied in the context of random graph 
theory |l|3j where it has been shown that connecting the vertices at random 
results in exponential degree distributions. The structure of a random graph is 
thus uniform, with most vertices having approximately the same degree and only 
few vertices of high degree. 

Empirical results indicate however that several real-world networks display 
instead power-law degree distributions PQ p{k) = cfc -7 , where c and 7 are con- 
stants. Because of the self-similarity properties of the power law, these networks 
are called scale- free. It is very remarkable that both the empirical and the analyt- 
ical work usually recovers scale-free networks with power law exponents ranging 
between 2 and 3. 

Since power laws decay much slower than exponentials, scale-free topologies 
are not uniform and display many highly connected vertices, called hubs. Hubs 
work as shortcuts between apparently distant environments and thus have a 
crucial role in many phenomena, such as in the spreading of infective agents 
through a network [3] or resilience of the network to external attacks [5] . 

Scale-free topologies are normally interpreted as a consequence of network 
self-organization, resulting in properties of the system as a whole which (in con- 
trast with classical physics) cannot be inferred only from properties of its parts 
(that is from the interactions between the individual components considered in 
isolation). 

Thus one of the major problems in complex network theory has been to 
identify the range of mechanisms by which a network can self-organise into a 
scale- free state [6]. 

It has been shown analytically |l|7|8|9j that if a network evolves by addition 
of new vertices at a constant rate and if the newly introduced vertices connect 
preferentially (and linearly [7j) to highly connected vertices, then the network 
displays a scale-free topology. 

Preferential attachment, however, is not always a natural hypothesis [519110] . 
For example, in many situations it is not realistic to assume that a node has the 
complete information about the degree distribution it would need in order to 
know where to attach preferentially. From the other hand this information is just 
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what is needed to normalise properly [34] the attachment probabilities as defined 
in most of the analytic models introduced in literature so far |l|6j ). Alternative 
mechanisms generating scale-free topologies have then been introduced, as for 
example the static model [T5], the varying fitness model |10I11I12I13I14| and 
random walk models incorporating a copying mechanism [34j . 

In addition to random networks, some interesting examples of determinis- 
tic networks have also been recently studied. A remarkable example are the 
Apollonian networks discussed in |16|17|18|19j , inspired by the ancient problem 
of finding a space filling packing of spheres. Apollonian networks turn out to 
be simultaneously scale-free, small world, Euclidean, space-filling and matching 
graphs. In this case the scale-free topology, however, turns out to be implicitly 
related to a preferential attachment mechanism [18] . 

It is important to recall that deterministic networks are, at least in principle, 
completely controllable and can thus be very useful in the design of artificial 
networks such as communication or economic networks, where (depending on 
the particular problem under consideration) the topology will usually have to 
satisfy several spacial and temporal constraints. 

The integer networks introduced in [20] are another interesting example of 
deterministic networks where vertices represent positive integers and edges are 
drawn whenever there is a divisibility relation between them. These networks are 
directed and their topology crucially depends on the particular set of vertices 
under consideration. For instance, when the vertices range in the set of the prime 
numbers an infinite star network centered in 1 is obtained. Considering instead 
the set of composite numbers, the degree distribution is the sum of an in-degree 
distribution and an out-degree distribution, and numerical simulations indicate 
that the latter is well approximated by a power-law with exponent 2. 

Thus our attention is drawn to abstract networks whose vertices are math- 
ematical objects and whose edges symbolise mathematical relations. Examples 
of this kind obviously abound. A straightforward example is the set of the in- 
finitely differentiable functions defined in the unit interval. We identify functions 
(through an equivalence relation) when they differ by a constant and naturally 
define a network structure on the set of equivalence classes by connecting vertex 
/ (for the sake of notational economy here we denote the equivalence classes by 
their representants) to vertex g if g is the derivative of /. Each non-constant ver- 
tex / has thus a descendent (its derivative /') and every vertex has an antecedent 
(its integral F(x) = f Q f(t)dt). We thus have defined a directed network con- 
sisting of infinite disconnected components. The degree of each vertex can only 
take the values 1 or 2 and the components can be described as infinite branches, 
finite branches (such is the case when one of the vertices of the branch is a poly- 
nomial), cycles (such as the one generated by sinx) and a loop (generated by 
exp(x)). 

Another way to associate a network to mathematical objects is by associating 
networks to equations. Consider any equation F(xi, X2, x n ) = 0, where the 
variables x±,X2, --^Xn range in a given set S. We can then naturally associate 
a network to the equation by representing the points of S as vertices and by 



4 



C. Bedogne', A. P. Masucci, G.J. Rodgers 



forming an n-clique (that is a completely connected subgraph) each time that 
the n elements of S satisfy the equation. We will then say that the network thus 
defined is generated by the equation F(xi,X2, ■■■,x n ) = 0. 

Cliques have been widely studied in the last few years [21j and often provide 
interesting insights into network structure and function. For instance cliques 
naturally represent clusters, communities and groups in social networks. Cliques 
are also very important in theoretical biology when modelling protein-protein 
interaction networks [22| and gene regulatory networks |23j . 

In this paper we will study the degree distribution of networks generated by 
Diophantinc equations, that is polynomial equations whose variables range in 
the set of integer numbers (originally introduced by Diophantus of Alexandria), 
attempting thus to establish a connection between age old unsolved problems in 
number theory and a very actual and practical area of Science. 

In particular, in Model A we study the finite networks generated by the 
Pythagorean equation x 2 +y 2 = z 2 when the variables range in the set {1, 2, N} 
recovering (at least up to N = 10 5 ) a degree distribution that is well approxi- 
mated by a power law with exponential cut-off. An ageing analysis reveals then 
that the hubs of the network, which are the nodes with the highest degree, are 
not the oldest vertices of the network as we would expect for a network grown via 
preferential attachment. In the Pythagorean case hubs form instead in vertices 
whose age is between old and middle age. After a while a freezing effect takes 
place and the hubs stop attracting new connections. Younger vertices then take 
over and become the new hubs in the network. This situation looks more akin to 
real social networks than to the rich get richer networks. This idea has already 
been suggested in a recent work in which rules are added to a stochastic 
network. The basic suggestion is that people in social networks agree and form 
hubs in response to an exact need, a deterministic rule, and not merely to an 
unbiased attraction to important people. The clustering coefficient analysis also 
displays interesting properties. In particular in Pythagorean networks c(k) has 
a power law behaviour that can be predicted analytically. 

In Model B we introduce the finite networks generated by the Diophantinc 
equation x 2 + y 2 = z in the set {1,2,..., N}. We analyse the degree distribution 
and clustering coefficient comparing the results to those obtained in Model A. 
In particular we find that, although the network's topology turns out to be 
very different than the one obtained starting from the Pythagorean equation, a 
power-law can still be recovered for several decades of values of k. 

2 Model A 

We consider the network generated by Diophantine equations of the form 

x n + y n = z n ^ 

where the variables x, y, z and the exponent n are positive integers. We study 
how the network generated by Eq. 1 evolves by introducing a new vertex at each 
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time-step, following the natural order of the positive integers. Thus at time- 
step t = 1 the vertex 1 is introduced and no edge forms. At each subsequent 
time step another integer is introduced and 3-cliques form whenever triplets of 
integers satisfy Eq.l. 

We then study the finite networks obtained by limiting the variables to range 
within sets of the kind {1,2,..., N}. The network will thus consist of N vertices 
(associated with the first N positive integers) and a number of edges belonging 
to 3-cliques. Notice first of all that at each time step the degree of every vertex is 
expressed by an even number. This fact alone has some interesting consequences 
in terms of network topology. For instance, it follows from a well-known theorem 
of graph theory that each connected component of the finite networks we just 
defined has an Eulcrian circuit, that is a circuit joining all edges and traversing 
each edge only once. Notice also that, according to Fermat's Last Theorem, 
proven by Wiles in 1994, Eq.l has no solutions for n > 3 and the corresponding 
networks are thus empty. We therefore need to consider only the cases n = 1 
and n = 2. 

2.1 The linear case 

We consider first the linear case 

x + y = z (2) 

with z < N. 

In this case it is a simple matter to derive the degree distribution analytically. 
Since multiple links are prevalent we performed a weighted analysis. We recall 
that the weight of the links is defined as the number of times two vertices are 
connected and the strength of a vertex as the total number of vertices it is 
connected to, having counted the repeated links. Since we allow loops, derived 
by a solution of EqJ5] of the type x + x = 2x, every vertex is connected to 
all other vertices if 2x < N, and to all other vertices but itself if 2x > N, so 
that k(x) = N - 1 if x < f and k{x) = N - 2 if x > f . If N is even we 
have P{k) = 2 " N ~_ 2 1) 5{k - N + 1) + 2 (N-i) S ( k -N + 2).HN is odd, P{k) = 
\5{k - N + 1) + \5{k - N + 2). 

To calculate the strength of the vertices we have to consider that, if n < 
^, the vertex n will assume 2(N — 1 — n) + 1 links (note that in this case 
the clique (n, n, 2n) will give 3 new links to vertex n ), forming the following 
N — 1 — n cliques: (n, 1, n + 1), (n, 2,n + 2)..., (n, TV — 1 — n, N — 1), plus n — 1 
links with the following cliques: (1, n — 1, n), (2, n — 2,n), (-|, ^, n) if n is even 
and (l,n - l,n),(2,7i- 2, n), (*±i - ^ n ^ L ,n) if n is odd. If n > f the 
vertex n will assume 2{N — 1—n) links forming the following N — 1 — n cliques: 
(n, 1, ri+ 1), (n, 2, ri+2)..., (n, N— 1 — n, N — 1) plus n — 1 links with the following 
cliques: (1, 7i— 1, n), (2, n — 2, ri), •••,(§, §, ri) if n is even and (1, n— 1, ri), (2, n — 
2, n), (^iii — 1, ^ji, n) if n is odd. It follows that the strength of a vertex n is 
given by s(n) = 2N - n - 2 if n < § and by s(n) = 2N - n - 3 if n > f . Thus 
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linearity of Eq.2 results in a uniform probability distribution for the strength: 
P(s) = j? if N - 3 < s < 2N - 3 and P(s) = otherwise. 

Notice at last that the degree distribution, in the thermodynamic limit, be- 
comes trivial as each vertex's degree diverges when TV — ► oo. 

2.2 The Pythagorean case 

We now study the network generated by the Pythagorean equation 

2,2 2 /r>\ 

x + y = z (3) 

Following standard mathematical terminology, we will say that (x, y, z) is a 
Pythagorean triplet whenever x,y and z satisfy Eq.3. The numbers x, y and z 
are usually called the "legs" of the triplet. When all legs are positive we say that 
(x, y, z) is a positive Pythagorean triplet. A Pythagorean triplet is called "primi- 
tive" whenever x, y and z are relatively prime. It has been known since antiquity 
that every Pythagorean triplet is either of the form (A(a 2 — b 2 ), 2\ab, A(a 2 + b 2 ) 
or of the form (2Xab, A(a 2 — b 2 ),X(a 2 + b 2 )). It is also well known that every 
positive integer x > 2 belongs to at least one Pythagorean triplet [27] and that 
there are infinite primitive Pythagorean triplets but, despite the apparent sim- 
plicity of the problem, the topological distribution of the set of all Pythagorean 
triplets is still largely unknown [24I25I26I27I28I29] . 

It is known for instance that the number A(N) of Pythagorean triplets with 
legs smaller or equal to a given positive integer N can be expressed by A(N) = 
^NlogN + BN + E(N) where B is an explicitly given constant and E(N) is a 
remainder function. It is also known that E(N) = 0(N?) and there is a large 
literature aimed at giving a better estimate for the remainder [24 25 26 27 2 8129] . 
It was proved in [S3] that E(N) = 0(N%exp{-c(logN)§ (loglogN)i) with some 
c > 0. Notice however that, in order to investigate the Pythagorean network 
just introduced analytically and derive its degree distribution we would clearly 
need to know much more, that is the exact number P(x) of Pythagorean triplets 
of which a: is a member. As far as we are aware this is still an open problem 
in number theory, probably deeply connected with the factorization problem. 
Indeed it is easy to conjecture that the properties of P{x) don't depend only on 
the algebraic form of Eq.3 but also on the structure of the number x itself in 
terms of its factorization as product of prime numbers. For instance it is known 
that if uu(x) is the number of distinct prime factors of an integer number x > 1, 
then the number P(x) of Pythagorean triplets of which x is a member is given 
by P(x) = 2 W ^ if all prime factors of x satisfy p = lmod4,P{x) = 2 W ^- 1 
if x is odd and some prime factor of x satisfies p = 3mod4, P(x) = 2 W ^ X ^ 1 if 
n = OmodA and P{x) = otherwise. The problem of determining the distribution 
of Pythagorean triplets is thus reducible, at least within the set of primitive 
triples, to the problem of counting the prime factors of each leg, which strongly 
suggests that there is little hope to derive the degree distribution analytically at 
present. The model, however, can be described numerically, at least in the finite 
case. We will then consider the networks generated by allowing the variables of 
Eq.3 to range only in the first N positive integers. 
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Numerical analysis Wc performed numerical simulations of the network growth 
with the deterministic attachment rule given in Eq.3. We show the results for 
the degree distribution in the left panel of FigQ]for z < 10 , 10 6 . What is first 
evident from the log- log plot is that, changing the scale of the network size, the 
degree distribution preserve its shape. Then, with increasing values of z, a wide 
region emerges where the probability distribution looks like to follow a power 
law. On the right panel of FigQ]we show the empirical distribution for z < 10 6 
after a logarithmic binning to understand better the function behind the data. 
Wc find that a power law with an exponential cut-off, 

P(k) oce-Afc" 1 ' 7 , (4) 

approximates very well the empirical data. 

Considering that complex growing networks usually display pure power law 
distributions only in the thermodynamic limit [24j . we find the finite results 
expressed graphically in Fig.l particularly remarkable. 




Fig. 1. In the left panel we show the degree distribution for the Pythagorean 
network for z < 10 , 10 6 . In the right panel we perform a logarithmic binning on 
the degree distribution for z < 10 6 to illustrate the power law behaviour. The 
fitting curve is the power law with exponential cut-off in Eqf?] 



To better understand how those structures form, we performed an ageing 
analysis of the network. In Figj2] we show the degree of the vertices versus their 
labels for different values of z. Small labels indicate old vertices. Peak values of 
the degree k represent the major hubs of the network. For stochastic preferential 
growing networks [30 , the high degree vertices are always the old vertices of the 
network. In the Pythagorean network this result does not hold. From Fig[3] we 
can see that after a while the old vertices' degree freezes and hubs form between 
young and middle-aged vertices. 
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Fig. 2. Representation of the aging of the hubs for the Pythagorean network for 
z < 1000, 2000 on the left and for z < 3000, 10000 on the right, k is the degree of 
the vertices, while n is their label. Small values of n correspond to old vertices. 



Since the Pythagorean network is made up of cliques it is interesting to look 
at its clustering coefficient. The clustering coefficient jl] is defined as 

where fcj is the degree of vertex i and e, is the number of nearest neighbors 
of vertex i that are connected to each other. We show in Figj3] our measures 
of c(fc) on a Pythagorean network for different values of z. The regularity of 
the distribution is striking and due to the high regularity of the network. All 
vertices with k = 2 form a triangle with their nearest neighbors, so that they 
have c(k) = 1. The higher a vertex degree, the more improbable it is that its 
nearest neighbors are fully connected. Since for the inner geometry of the network 
e i > % ' ^ f°ll° ws that c(fcj) > and, obviously, c(ki) < 1. All the values of 
c{kj) lie on well-determined levels defined by 

*» = ^)' (6 > 

where n = 0, 1, 2, rij and rii = k ^ k ^~ 2 > . Xo demonstrate Eqj6] is a simple 
graphical exercise. Eqj6] is compared to the numerical simulation in Fig[3j For 
each ki, the different values of c(fcj) are equally spaced at a distance of A(c(ki)) = 

2 

k<(k«-i) ' 

Notice at last that, in contrast with the linear case, the network's behaviour 
in the thermodynamic limit is now highly not trivial and probably cannot be 
described in detail before some open questions in number theory will be even- 
tually answered. The numerical work we performed may however suggest some 
interesting conjectures. We think that the most relevant feature of our model 
is perhaps the ageing analysis showing that, when the network grows, the hubs 
are constantly rejuvenated, in sharp contrast with the behaviour of scale-free 
networks generated through the preferential attachment mechanism. 
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Fig. 3. Clustering coefficient for the Phytagorean network measured in simula- 
tions with z < 1000, 10000. The lines represent EqJB]for different values of n. In 
particular the greatest value of n allowed for z < 10000 is n = 23. 



We now study the network generated by the deterministic rule expressed by the 



This network is thus deeply connected with the classical Euler and Fermat Prob- 
lem of finding all integers that can be expressed as the sum of two squares. In 
a sense this case is intermediate between the linear and the Pythagorean case. 
Notice, first of all, that in the thermodynamic limit EqfT] generates a trivial net- 
work where each vertex's degree goes to infinity. Indeed given any positive integer 
x, for any arbitrary integer y (when N — > oo) wc can always find an integer 
z satisfying EqfT] Thus we only need to consider the finite networks specified 
by the inequality z < N. The network attachment mechanism, as in Model A, 
will depend on the generating equation in the sense that every triple of integer 
numbers obeying EqfT] forms a clique in the network. 

The resulting network is very well connected since every couple of numbers at 
a certain time step (when N is sufficiently large) forms a clique satisfying Eq\7\ 
Nevertheless interesting topological properties emerge. In the left panel of FigfJ] 
we show the degree distribution for this network. The shape of the distribution 
is well approximated by a power law with exponent —4.5 for several decades of 
values of k. Very interestingly, after falling the distribution begins to rise again. 
On the right panel of the same figure we show the degree of the vertices against 
their label number n. This number can be considered as the birth time of the 
vertex so that wc can observe the evolution of the hubs of the network with 
ageing. The degree is nearly constant and very high for small values of n and it 
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equation 



x 2 +y 2 = z 



(7) 
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drops rapidly to small values for n ~ \/N (this comes from Eq.7, when z ~ N 
then x, y ~ 

Notice that the power law exponent is considerably lower (at least from 
the point of view of network theory where power law exponents usually range 
between —3 and —2) than the one we recovered in the Pythagorean case. We 
interpret this weakening as the result of the upper thresholds we imposed to the 
network which, in this case, have a much stronger effect on the degree distribution 
than in the Pythagorean case. 

In Figf5]we show the clustering analysis for the network. On the left panel 
we show the clustering coefficient c(k), defined in Eqj5l against the degree, while 
on the right panel we show the average clustering coefficient < c(k) > against 
the degree. We can observe that it is very high for small values of k, where the 
degree distribution is a power law, then it is nearly constant for several decades 
and finally it falls rapidly to zero for high values of the degree. This behaviour 
has to do with the clique structure of the network. c(k) = 1 as soon as isolated 
cliques form, that is for k = 2, then its value falls as new connections are added 
and the degree increases. 




10 100 1 10 100 1000 10000 100000 

k n 



Fig. 4. Left panel: degree distribution for the deterministic network generated 
by the rule in Eq{7|for z < 1000, 10000, 100000. Right panel: representation of 
the aging of the hubs for the same network for z < 10000, 100000. k is the degree 
of the vertices, while n is their label. Small values of n correspond to old vertices. 



Conclusions 

We have introduced the networks generated by Diophantine equations focus- 
ing on two classical problems of number theory widely studied since antiquity. 
We investigated first the Pythagorean problem of finding all triples of positive 
integers that can be represented geometrically as legs of a right-wristed triangle. 
We studied then the (closely related) problem of determining the set of positive 
integers that can be represented as the sum of two squares. These two ques- 
tions, despite their apparent simplicity, have not yet been answered analytically 
in their full generality. 
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Fig. 5. Left panel: clustering coefficient c(k) against the degree k for the de- 
terministic network generated by the rule in Eq|7]for z < 10000. Right panel: 
average clustering coefficient c(k) against the degree k for the deterministic net- 
work generated by the rule in EqJ7]for z < 10000. 

We noticed that these classical problems can be naturally translated into the 
language of network theory and studied them numerically in their finite versions 
(that is by not allowing integers to be greater than a fixed upper threshold). 
It is then easy to see that the equations considered generate clique structured 
networks, in an interesting analogy with many relevant real world networks such 
as social [31] and biological networks [22|23j . 

In the Pythagorean case, where the network evolves through the deterministic 
attachment mechanism specified by Eq.3, we recovered a degree distribution 
consistent with a power law with exponential cut-off, in accordance with the 
behaviour displayed by many real world random networks PQ . 

We then performed an ageing analysis and showed that, in sharp contrast 
with stochastic preferential growing networks, hubs form between young and 
middle aged vertices (the degree of each vertex, after a while, freezes). 

The power-law, in this case, is clearly not a by-product of any preferential 
attachment or fitness mechanism, but a consequence of the algebraic structure 
of Eq.3 itself. This shows once more that preferential attachment is a sufficient 
but not a necessary condition to the generation of scale-free topologies. As far 
as we are aware the network generated by Eq.3 is the first genuine example in 
literature of a deterministic scale-free network (at least in some approximation) 
which does not grow preferentially. 

Power laws are normally regarded as a sign of complexity and we found it 
particularly stimulating to detect them not only in natural phenomena but also 
at the roots of number theory itself, indicating a fascinating connection between 
complex systems and pure mathematics. We then studied the second classical 
problem mentioned above studying the network generated by Eq[7]and recovered 
a degree distribution approaching a power-law with a relatively high exponent 
7 ~ 4.5 for several decades of values of k. 

We believe that the study of networks generated by equations will turn out 
to represent an important field of interdisciplinary research and we hope that 
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our discussion will be useful to researchers in both complex network theory and 
pure mathematics as a first step to recover more rigorous and general results. 
In particular, in future work it will be stimulating to investigate the relation 
between an equation's algebraic properties (or analytic properties if, for instance, 
differential equations are considered) and the associated network's topology. 

Diophantine networks may find interesting theoretical and technological ap- 
plications. First and foremost they can be used to design deterministic toy models 
of complex systems, allowing to find practical ways to build networks system- 
atically without having to deal with different degrees of stochasticity in their 
architecture (in this sense they could play a role similar to the one played by 
the Euclidean grid in other traditional contexts). Diophantine networks may 
also allow to overcome some limitations intrinsic in the preferential attachment 
method, namely the fact that the hubs are usually the oldest vertices and never 
"die", leaving little room for network rejuvenation. Here instead we have a dy- 
namic situation where hubs are dominant, during network growth, only for a 
limited period of time. There may be situations where this feature is realistic 
(for instance in models of technological developments). Another interesting ap- 
plication of the theory may be found in the study of the behaviour of dynamically 
evolving agents. 
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